Grain refinement and partitioning of impurities in the grain 
boundaries of a colloidal polycrystal 
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We study the crystallization of a colloidal model system in presence of secondary nanoparticles acting as impurities. Using con- 
focal microscopy, we show that the nanoparticles segregate in the grain boundaries of the colloidal polycrystal. We demonstrate 
that the texture of the polycrystal can be tuned by varying independently the nanoparticle volume fraction and the crystallization 
rate, and quantify our findings using standard models for the nucleation and growth of crystalline materials. Remarkably, we find 
that the efficiency of the segregation of the nanoparticles in the grain-boundaries is determined solely by the typical size of the 
crystalline grains. 
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1 Introduction 

Impurities affect drastically crystal nucleation, growth and 
texture and understanding the role they play during the solid- 
ification processes is of great interest in material science and 
engineering l , pharmaceutical industry 2 , mineralogy 3 , protein 
crystallization 4,5 and life science 6 . 

Moreover, the control of the texture by addition of solutes 
is important for tailoring the mechanical properties of crys- 
tals^. 

Relevant progress in understanding crystallization of mono 
and bidisperse systems has been achieved using colloids, that 
are often regarded as an analog to atoms on larger length and 
time scales 10 , enabling real time studies by direct visualiza- 
tion of processes of importance and interest in condensed mat- 
ter physic s u i 12 . Crystallization of nano- and micro-spheres 
have been widely investigated numerically 13 and experimen- 
tally 14_16 , including in bi-component systems. Binary mix- 
tures of colloids where the volume fraction of the two species 
is not too unbalanced (up to a factor of 7, but most typically 
[by a factor of 2) exhibit a great variety of crystalline super- 
lattice structures similar to atomic systems 17_2 °. By contrast, 
the addition of a very small amounts of a dopant in a solidi- 
fying matrix does not change the crystalline unit cell, but can 
modify the nucleation—"—, and induce crystal growth frustra- 
tio n 24 ' 25 local defect formation 26 and local fractionation 27 , 
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depending on the volume fraction and size of the dopants. 

In this article we present a novel strategy for the control of 
the texture of colloidal polycrystals, based on the addition of 
small amounts of nanoparticles (NPs) as dopants to a solidi- 
fying matrix. We achieve the segregation of NPs in a network 
of thin sheets that is the colloidal counterpart of the texture 
of grain boundaries formed during the solidification of crys - 
talline metallic alloys 28 and other molecular materials doped 
with impurities 2 -. By varying the NP volume fraction and 
the rate at which the background matrix is solidified, both the 
typical mesh size of the dopant network and the efficiency of 
the partitioning of the NPs between the matrix and the grain 
boundaries are tuned. In the framework of the classical nu- 
cleation theory, we show quantitatively for the first time how 
both the nucleation and the growth of colloidal crystallites are 
influenced by the NP volume fraction and the crystallization 
rate, determining the final grain size. 



2 Materials and methods 

The colloidal polycrystal is composed of an aqueous solution 
of Pluronic F108, a commercial PEO-PPO-PEO triblock co- 
polymer (Serva Electrophoresis GmbH), where PEO and PPO 
denote polyethylene oxide and polypropylene oxide, respec- 
tively. The co-polymer (concentration of 34% w/w) is fully 
dissolved at T ~ °C. Upon increasing T, it self-assembles 
into spherical micelles of diameter d ~ 22 nm, due to the 
increased hydrophobicity of the PPO block 30 . The volume 
fraction, <j>, of the micelles increases with T, until crystalliza- 
tion occurs due to micelle crowding. Differential scanning 
calorimetry and rheology show that crystallization occurs be- 



tween 15 °C and 16 °C, depending on the heating rate—. We 
emphasize that, unlike conventional colloidal systems where 
controlling precisely and changing in situ the volume fraction 
is a difficult task, our system allows crystallization to be in- 
duced at the desired rate simply by varying T, in analogy with 
molecular materials. 

In order to analyze in a quantitative way the crystallization 
process, as it will be discussed in Sec. [3] it is convenient to 
map our micellar system onto a hard sphere suspension, by 
determining the relation between T and (p. The soundness of 
this approach is supported by previous work 32 showing that 
the static structure factor of a similar micellar system could 
be mapped onto that of hard spheres using an affine relation 
between T and (p. In our case, we determine <j>(T) by com- 
paring the temperature dependence of the sample viscosity in 
the micellar fluid phase to the volume fraction dependence 
of the viscosity of hard-spheres suspensions 33 . The mapping 
thus determined is <j) — a x (T — To) with Tq — 4.2 °C and 
a ~ 0.045 °C~ 131 . With this mapping, rheology measure- 
ments in the regime where the sample is solidified very slowly 
so as to avoid any significant undercooling show that the on- 
set of crystallization occurs at = C tyst ~ 0.5136, a value 
compatible with that expected for slightly polydisperse hard 
spheres. 

We dope the micellar crystal with small amounts (at most 
1% v/v) of green-yellow fluorescent carboxylated polystyrene 
nanoparticles purchased from Invitrogen, with diameter a = 
20,36 and 100 nm. Neutron scattering measurements show 
that the microscopic crystalline structure (face-centered cubic, 
with lattice parameter 30 nm) is preserved upon addition of 
up to 2% v/v of NPs- 3 -^, twice than the highest concentration 
used here. For confocal microscopy imaging, samples are in- 
troduced in chambers of thickness 250/xm and solidified by 
raising the temperature from 3°C to 23°C at a controlled rate, 
t , and visualized at least 20/xm from the walls. 

3 Results and Discussions 

We show in Figs. [T]\-C representative images of colloidal 
samples doped with NPs with diameter a = 36 nm. All sam- 
ples are solidified at the same rate, t = 7 x 10~ 3 °C min -1 , 
and the average NP volume fraction cq is varied between 
0.05% and 1%. Images exhibit fluorescent (clear) connected 
thin lines that separate featureless (dark) zones. Confocal 
scans through the material show that the clear lines are in fact 
2-D sections of a network of grain boundaries enriched in NPs 
that delimitate grains with different crystalline orientations. 
A 3-D reconstruction of sample with co=0.5% is reported in 
Fig.QJ), where different false colors indicate different grains. 
FiguresQ~|\-C clearly show that the grain size decreases by in- 
creasing the amount of NPs. To quantify the influence of NPs, 
we calculate the cross section area of the grains, S. The fre- 



quency distribution of S for various cq is reported in Fig. [TJi 
[about 800 cross sections, taken in distinct positions on both 
sides of the chamber are measured for each experimental con- 
dition]. The peak of the distribution moves to smaller values 
with increasing co, demonstrating grain refinement at larger 
impurity content. 




Fig. 1 Color online. (A-C) Confocal microscopy images of three 
samples prepared at different impurities (NP) volume fraction cq. 
The clear lines are sections in the observation plane of the grain 
boundaries, where the NPs accumulate. The temperature ramp is 
f — 7 x 10~ 3 °C min -1 . (D) 3-D reconstruction of a sample with 
cq = 0.5%. The size of the reconstructed region is 35 x30x5 /xm 3 . 
(E) Frequency distribution of the cross section of the grains, for 
samples with various co prepared at a constant f = 
7x 10 °C min . Inset: same data as in the main plot but scaled 
by the average cross section S. Impurities are fluorescently labeled 
polystyrene NPs with diameter 36 nm. 

In analogy to atomic crystals and other materials 35 , we find 
that the rate at which temperature is changed during solidifica- 
tion has also a crucial influence on the final grain size, slower 
ramps leading to a coarser texture. This is demonstrated for 
samples with a fixed NP volume fraction (cq = 0.5%) and 
prepared with various temperature ramps, as in the images 
in Figs. |2j\-C and in figure |2jD where the size distributions 
of the cross section area of the grains, spanning two decades, 
are reported. In all cases, the distributions are relatively nar- 
row, their standard deviation normalized by the average being 
0.7 ± 0.2. Remarkably, we find that the shape of the distribu- 
tions is almost independent of the NPs volume fraction and of 
the temperature ramp, as demonstrated by the inset of Fig. [TJi 
(resp. Fig.|2jD), where curves for various cq (resp. f) are col- 



lapsed by normalizing the S with respect to its mean value, S. 
In the insets, the right tails of the distributions are approxi- 
mately straight lines in a semi-logarithmic plot, indicating ex- 
ponential tails. 
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Fig. 2 Color online. (A-C) Confocal microscopy images of three 
samples prepared at various t. The NP volume fraction is fixed at 
Co = 0.5%. (D) Frequency distribution of the cross section of the 
grains, for samples with Co = 0.5% and prepared with various t . 
Inset: same data as in the main plot but scaled by the average cross 
section 5. Impurities are fluorescently labeled polystyrene NPs with 
diameter 36 nm. 

Figures [5K and B display the evolution of the mean ra- 
dius of the grain, R = y S/lt, where S is the mean cross- 
section area, with impurity amount cq and temperature ramp 
f, respectively. The grain size is of the order of a few 
micrometers^, with a clear influence of both NP content and 
crystallization on the size. 
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Fig. 3 (A) Mean radius of the grains vs. volume fraction of 
impurities, at a fixed temperature ramp f= 7x 10 °C min and 
(B) vs. f, at a fixed impurity volume fraction co=0.5%. Symbols are 
experimental data and lines are numerical calculations, as explained 
in the text. The dashed line is the theoretical * dependence in the 
absence of impurities. In A and B the bars are the standard deviation 
of the size distributions. 



' Note that strictly speaking R is smaller than the actual average grain radius, 



To quantitatively analyze our experimental findings, we use 
classical models for the nucleation and growth of crystalline 
grains. We first recall how an average grain size can be com- 
puted, using the Johnson-Mehl-Avrami-Kolmogorov (JMAK) 
theory 36-40 which assumes that nucleation occurs randomly 
and homogeneously, and that the growth velocity of the crys- 
tallites, v g , does not depend on the extent of the crystallization 
process and is isotropic. Following JMAK, the extended vol- 
ume fraction of the sample that has crystallized at time t , X e (t), 
reads: 
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where /(t) is the nucleation rate per unit volume at time T. 
The extended volume fraction does not take into account that 
as crystallization proceeds less fluid phase is available for nu- 
cleating new crystallites. Moreover, it doesn't account for the 
fact that crystallites stop growing when they impinge into each 
other. As shown by Kolmogorov and Avrami, these effects can 
be accounted for by calculating the actual volume fraction that 
has crystallized, X a (t), according to : 



X a (t) = l-exp[-X e {t) 



(2) 



The expected grain radius R can be calculated from the 
grain number density at the end of the crystallization pro- 
cess 41 , yielding 
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where l a is the actual nucleation rate: 

4(t) = [1-X„(t)]/(t) 



(3) 



(4) 



Note that the nucleation rate, /, and growth rate, v g , are in 
general time-dependent quantities. This is indeed the case in 
our experiments, where crystallization occurs while the sam- 
ple is submitted to a continuous increase of temperature, lead- 
ing to a continuous increase of the micelle volume fraction. 
In the following, we aim at calculating R using Eqs. (Q]|4| in 
order to compare it to the experimental results of Fig. HJ^-B. 
To this end, explicit expressions are required for / and v g . We 
start by discussing the growth velocity and assume that v g is 
given by the Wilson-Frenkel law—: 



, = ^[l-exp(-|A-i|/fe B r)] 



(5) 



where k% is Boltzmann's constant, Ajj, is the difference in 
chemical potential between the solid and the liquid phases, 

since S is the grain cross section obtained by a random slicing, not the max- 
imal cross section. However, we expect the difference to be relatively small: 
for an isolated sphere, e.g., the average radius of a random cross section dif- 
fers from the actual sphere radius by less than 20%. 



D s is the self diffusion coefficient of the micelles in the liquid 
phase, and A is a typical distance over which diffusion occurs 
when a micelle in the fluid phase is incorporated in the grow- 
ing crystallite. 

To evaluate Ajj,, we take advantage of numerical and ex- 
perimental results on the crystallization of hard spheres (HS) 
suspension s 13 i 14 i 16 i 42 ; assuming that micelles in first approx- 
imation behave as HS and using the (j)(T) mapping dis- 
cussed in Sec. [2] We evaluate Ajj. from numerical simu- 
lations of HS colloids 13 , where data can be well fitted by 
AjU =A( l kBT(<j> C ry S t - 0) with <p crys , w 0.5 and A M sa 15. In the 
following, we take C ryst =0.5 136 as discussed in Sec. HI while 
Ap is treated as an adjustable parameter. The volume fraction 
dependence of D s for the micelles is estimated from experi- 
mental data on the structural relaxation time associated with 
self-diffusion in suspensions of HS colloids 43 , assuming that 
the diffusion coefficient is inversely proportional to the struc- 
tural relaxation time. We use the empirical formula D s =Dqx 

(1 - 2.50 + 1.360 2 )(1 - $;) + 144.570 2 (1 - ^) 2 - 5 ] , which 
reproduces very well the experimental data in the range < 

6 x 10~ 14 m 2 /s as the dif- 



(j) < 0.59, and use D 



k b t 



fusion coefficient of a micelle in the dilute regime, where 
Tjo = 0.32 Pa s is the solvent viscosity. The quantity A is ex- 
pected to be of the order of the micelle diameter, d, although 
significative differences between simulations and experiments 
have been reported—. We set A = A\d with Ax an adjustable 
parameter. 

We now turn to the nucleation rate I. Classical nucleation 
theory (CNT)ii predicts 



/ = rexp[-AG*/Jt B r] 



(6) 



where AG* is the nucleation barrier in the Gibbs free energy 
to be overcome in order to form a stable nucleus, and the ki- 



netic prefactor r = ( -. 



,^^|A M |) V^»' /3 . Here« c isthe 
number of micelles in the critical cluster, and pi = 6<j)/nd 3 
is the number density of the micelles in the liquid phase. Ac- 
cording to CNT, AG* is the maximum in the Gibbs free energy, 
which is the sum of a surface free energy term, AGs = 47Tr 2 yo, 
and a volume term, AGy = — ^-P ps\Afi\. Here r is the ra- 
dius of a nucleus, 70 is the surface free energy density of the 
fluid-crystal interface and p$ is the number density of micelles 
in the solid phase (p s — 6(j> c /nd 3 with <j) c = 0.5). One readily 
finds 



AG* 



167T Yo 



(7) 



3 p s |A^' 

For the surface free energy density, we use 70 = ksTAy/d 
with Ay an adjustable parameter expected to be of order 
1 n i 42 i 45 . Numerical simulations of HS— show that n c the 
number of colloids in the critical nucleus is of order 100 and 
depends only weakly on the volume fraction: accordingly, we 
take n c = 100 independent of (p. 



In order to model the effect of the addition of NP, we need 
to modify the CNT expressions. We assume that the most rele- 
vant effect is a modification of the surface free energy term due 
to the accumulation of NP at the interface between the crystal 
and the liquid phase (see fig. |5), since slight variations of 70 
are known to have a massive impact on nucleation rates—. We 
thus neglect any dependence of |A/i| on NP content, implying 
that both Vg and the volume term AGy are unchanged upon ad- 
dition of NP, and propose that the presence of the NPs at the in- 
terface leads to a reduction of the surface free energy density. 
We assume that this reduction is proportional to the amount of 
NPs at the interface, which in turn is proportional to n p , the 
amount of NPs initially available in the volume occupied by 
a crystallite. Hence, for a nucleus of radius r, the reduction 
of surface energy is proportional to n p = 8r 3 co/(7 3 , and to the 
cross-section of the NPs, no 2 /A. The surface term then be- 
comes AGs = 47Tr 2 yo "^" r 3 , with F p a proportionality 

constant. Note that the NP contribution, although physically 
due to a change of the surface term, formally can be incorpo- 
rated in the volume term of AG, because it scales with r 3 . In 
the presence of NPs, the free energy barrier is then obtained 
by a suitable modification of the denominator of Eq. ©: 
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Using Ajj, Ax, Ay and Fp as adjustable parameters, we fit 
the experimental R using the modified CNT approach out- 
lined above. Operationally, we take t — as the time at which 
(j) = ciyst and, for a given f and cq, we calculate numerically 
X e (t ) using Eq. (JTJ and v g (t) and 7(f) as given by Eqs. J5I6I8I ), 
where the latter quantities depend on time through D s and Ajj.. 
Equations (f2]|4|i are then used to obtain R. We show in Fig.[3}\- 
B the best fits as solid lines. While not perfect, the fits exhibit a 
reasonable quantitative agreement with the experimental data, 
especially given the order-of-magnitude variations of quanti- 
ties related to colloidal crystallization often observed in the 
literature^. The fitting parameters are: Ax = 0.14 ±0.05, 
A M = 17 ± 3, A r = 0.747 ± 0.20 and F P = 126 ± 24. The 
three fitting parameters that are independent of NP content 
and for which a comparison with previous works is possible 
are in good agreement with numerical results 13 , for which 
Ax = 0.2 - 0.5, A M « 15, Ay = 0.64. Furthermore, this simple 
model is able to capture the main experimental trends, in par- 
ticular the roll-off of R(f) at slow ramps, which is due to the 
presence of NPs (compare the continuous and dotted lines in 
Fig. |5J}) . The discrepancy at high cq observed in in Fig. fJK 
is presumably due to an over-evaluation of the effect of NPs 
on the surface tension, due to the fact that, as we shall show 
it in the following, as cq increases the fraction of NPs that 
are effectively expelled from a growing crystallite decreases. 
Although a qualitative relation between AG* and the typical 
crystallite size has been reported for metals 46 , our work is the 



first quantitative analysis for colloids of the grain size using 
nucleation and growth theories^. 

So far, we have shown that both the volume fraction of im- 
purities and the crystallization rate determine the polycrys- 
talline texture, thereby allowing for a fine control of the av- 
erage grain size over almost one decade. We now turn to the 
investigation of the degree of confinement of NPs once the 
solidification process is completed. To better appreciate the 
spatial distribution of the NPs, we show in Figs. |4j\-B con- 
focal images of samples doped with particles large enough 
(<r = 100 nm) to be individually resolved by confocal mi- 
croscopy, at least for moderate volume fractions. The dopant 
volume fraction is kept fixed at cq — 0.05%, while the effect 
of f is investigated. 
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Fig. 4 Color online. (A-B) Confocal images of samples prepared 
with ramp rates 7x 10~ 3 °C min -1 and 5x 10 -4 °C min~', 
respectively. (C) Ratio p of the volume fraction of NPs in the grain 
boundaries and inside the grains. The symbols are labeled by the 
particle diameter and their number density, p. (D) Same data as in 
(C) collapsed onto a master curve by plotting pasa function of the 
average grain cross section. 

As shown in Fig.[4}\, for polycrystals prepared using a rel- 
atively fast temperature ramp, an excess of NPs is observed 
in the grain boundaries, consistently with our previous ob- 
servations with smaller dopants. However, large numbers of 
individual NPs are clearly visible also in the interior of the 
grains. By contrast, NPs are almost completely expelled from 
the grains in samples prepared using a slower ramp (Fig. 5B). 
We quantify the efficiency of the partitioning by evaluating the 
ratio, p, between the NP volume fraction in the grain bound- 
aries and that in the interior of the grains, as obtained from 
fluorescence intensity measurements. Figure |4t shows that 
confinement in the grain boundaries becomes increasingly ef- 
ficient as the temperature ramp slows down, similarly to ob- 
servations in atomic systems 47 , and that, quite generally, p 
decreases when the average number density of impurities, p, 



increases. Thus, an increase of either f or p yields both a 
finer texture and a smaller degree of partitioning, suggesting 
that the degree of partitioning should correlate with the grain 
size. We note that the grain size, directly related to the number 
of successful nuclei, does not depend on the size of NPs but 
only on their quantity, signaling the absence of heterogeneous 
nucleation on the surface of NPs, as expected since the size of 
the NPs is comparable to the size of the micelles 21 . 
Remarkably, we find that all data for ramp rates spanning two 
decades, impurity number density spanning three decades, and 
for NP diameters ranging between 20 nm and 100 nm indeed 
collapse onto a master curve when plotting p as a function of 
the average grain cross section (Fig. HJ3). Interestingly, p is 
close to unity and rather insensitive to S until S sw 800 |im 2 , 
above which the degree of partitioning increases steeply. 

To understand the origin of the correlation between grain 
size and degree of particle segregation, we analyze the growth 
of the crystallites during solidification, a quantity not accessi- 
ble in metallic systems. This is illustrated in Fig.|5]and Sup- 
plementary movie S 1 . 

At low temperature, below the onset of crystallization, the 
images are homogeneously grey down to the smallest length 
scale that is accessible (1 pixel = 0.11 /im), revealing that 
the impurities are homogeneously distributed in the sample, 
which is still fully fluid. The nucleation and early stage growth 
can not be resolved in our confocal microscopy experiments, 
presumably because the partition of the fluorescent impurities 
is too weak (Fig. UK)- As solidification proceeds, however, 
the interface of the growing crystallites becomes increasingly 
visible, because the fluorescent NPs are expelled by the solid 
phase and accumulated in the fluid boundary layer (Fig. |5Ji- 
C). Eventually, the different grains that are growing touch each 
other and become facetted, resulting in the formation of a net- 
work of grain boundaries where the NPs are confined at a den- 
sity larger than in the bulk (Fig. EJ)). In Fig. [5ji we report the 
normalized impurity volume fraction profiles, obtained from 
the fluorescent intensity, along the growth direction of one 
grain from its center: the height of the peak of NP volume frac - 
tion increases during the propagation of the solid-fluid inter- 
face, indicating that the amount of impurities rejected from the 
crystallite and accumulated at the interface increases steadily, 
until the grain stops growing when impinging on an adjacent 
growing crystallite. Thus, the larger a grain can grow, the 
higher the volume fraction of the impurities eventually trapped 
in the grain boundary, explaining the correlation between the 
grain size and the degree of partitioning shown in Fig. |4jD. 

As a concluding remark, we note that we have successfully 
reproduced our experiments with a variety of NP of different 
kinds, thereby demonstrating that the system described here 
provides a powerful, robust and convenient way to study the 
crystallization of bi-disperse colloidal systems and to segre- 
gate NPs in thin sheets within a 3-D soft matrix. The colloidal 
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Fig. 5 Color online. (A-D) Confocal microscopy images showing 
the progressive accumulation of impurities at the fluid-solid 
interface during solidification. Images are taken at different 
temperatures, as specified by the labels, during solidification at a 
rate 7=0.03 °C min . The impurities are 36 nm fluorescent NPs, at 
a concentration of co=0.05%. The white rectangle shows the region 
of the growing crystallites used for the analysis of the temporal 
evolution of the concentration profile shown below. (E) Time 
dependence of the normalized NPs volume fraction obtained from 
fluorescence intensity profiles calculated across a solid-fluid 
interface. The doted vertical lines indicate the position of the 
solid-liquid interface as determined by the peak of the fluorescence 
profiles. 



composite material investigated here, thanks to its softness and 
optical transparency, offers a unique opportunity to investigate 
in real time and by direct visualization the formation of a poly- 
crystalline texture and its response to an external perturbation. 
In particular, we are currently exploring the dynamics of the 
grain boundaries under a mechanical load. 
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